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Entanglement entropy in d + 1 SU(N) gauge theory 
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We consider the entanglement entropy for a sub-system in d + 1 dimensional SU(N) lattice gauge 
theory. The 1 + 1 gauge theory is treated exactly and shows trivial behavior. Gauge theories in 
higher dimensions are treated within Migdal-Kadanoff approximation. We consider the gauge theory 
in the confinement phase. We demonstrate the existence of a non- analytical change from the short 
distance to long distance form in the entanglement entropy in such systems (d > 2) reminiscent 
of phase transition. The transition is manifested in nontrivial change in the RG flow of character 
expansion coefficients defining the partition function. 
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I. INTRODUCTION 

Interest in the study of entanglement entropy has a 
relatively long history in quantum field theory. The 
early motivation was due to the connection with black 
hole physics. General properties of the entanglement en- 
tropy, such as its dependence only on the surface, were 
demonstrated for system of oscillators and massless non- 
interacting scalar field theory [l[ . Later the entanglement 
entropy was studied in gravity duals of confining large N c 
gauge theories [2| using the AdS/CFT approach of [3j. In 
this work the d dimensional space was divided into two 
complementary regions A and A by two imaginary d — 1 
dimensional hyper-surfaces placed distance I apart along 
one of the space directions 



A 
A 



td-1 



X JLj, 



(1) 



where I; is a line segment of length I . The authors studied 
the entanglement entropy as a function of / and found 
that it exhibits a non-analytical chan ge i n behavior at 
I = I* reminiscent of a phase transition[22[. 

In the present work we aim to prove that this is a gen- 
eral scenario for SU(N) gauge theories (at arbitrary N) 
at temperatures corresponding to the confinement phase. 
We consider ad+l dimensional gauge theory at finite 
temperature T. The zero temperature system is recov- 
ered as the limit T — > 0. We consider the same geometry 
of the entangled region as in Q , see eq. |T]) . 

If we are given the density matrix for such a system, 
we can integrate out all degrees of freedom associated 
with region A. The resulting density matrix can be used 
to construct the entanglement entropy 



Pa — Tr^p, S A = -TiaPa log pa, 



(2) 



which is the entropy as seen by an observer with no access 
to the degrees of freedom in A. 



We will use a method of gluing replicas of the system 
under consideration into a multi-sheet Riemann surface, 
which was used in an extensive treatment of 2d CFT in 
[3, [1[ . We consider n replicas of a system after the trace 
over region A has been taken (for this the boundaries in 
time direction of this region were identified) [23[ . Each of 
these replicas is glued to another along the boundary of 
region A normal to the time direction. The first replica's 
upper boundary (coordinate t = 1/T) is glued to the 
lower boundary (coordinate t — 0) of the second replica 
and so on. The upper boundary of the last n-th replica 
is glued to the lowest boundary of the first replica, thus 
closing the system. For illustration of such gluing in 1+1 
and 2+1 dimensional theories see Figs. [T]and[21 One can 
observe that in such a system 



Trp A = 



Zn(A) 



(3) 



where Z n is the partition function of the glued system 
and Z is the standard partition function of the original 
system (Z = Z\). 

This approach allows one to construct the entangle- 
ment entropy 



S A 



lim ^-Trp n A 
n-»i dn ' A 



lim AMI) (4) 

dn Z n w 



II. SU(N) GAUGE THEORY IN d+1 
DIMENSIONS 

The partition function for SU(N) lattice gauge theory 



is 



Z = 



(5) 
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where the action is S p = S(U P ) = -f3/(2N)TrU p + h.c, 
j3 = 2N/g 2 is the lattice inverse coupling, and the pla- 
quette variable is the ordered product of gauge fields 
which live on the links constituting the plaquette U p = 



2 



ILedp The S au S e invariant action is a class function 
and therefore it can be expanded in group characters 

e -S p = ^ F r d rX r(Up) = F 1 + ^ CrdrXAUp) 
r \ r=iO 



where the first sum runs over all irreducible representa- 
tions, while the second sum excludes the trivial r = 
representation. For general SU (N) group r is a set 
of indices; d r is the dimension of the representation, 
c, r = F r /Fo < 1 and F r are the coefficients of expan- 
sion 



F = I dUe- s ^jxl{U). 



A. d — 1 gauge theory 



(7) 



R 



1/T 



FIG. 1: Z n for 1 + 1 dimensional gauge theory. 

The 2-dimensional SU(N) gauge theory is exactly solv- 
able, see @ for an overview and large N treatment of zero 
temperature U(N) gauge theory. It is possible to treat 
the zero temperature gauge theory in this study in anal- 
ogous fashion, however we would like to consider a more 
general case of a gauge theory at finite temperature T and 
therefore adopt a different approach. Finite temperature 
gauge theory in 2 dimensions normally is formulated on 
a R x §i surface periodic in time direction with period 
1/T. For practical reasons we consider a finite system in 
space direction. The corresponding discretized theory is 
formulated on a N r x N t lattice, with space-time cut-off 
a and aN t = 1/T and aN r = R. 

At this point it is instructive to consider how integra- 
tion on the surface is performed. Consider an elementary 
part of a surface bounded by a single loop 1241 . Its con- 
tribution to the resulting partition function[25j is 



f({a};dA) = l 



y^.diaiXijdA), 



(8) 



where OA is a product of link variables along the surface 
perimeter and the function is defined when all coef- 
ficients are specified. The expression for a junction of 
two surface elements A and B with a common boundary 
An B is 



(6) f({c};d(AUB)) 



d(AnB)f({a};dA)f({b};dB) 



l + J2di(kXi(d{AUB)), 



(9) 



The integration over the common boundary U — A n B 
is performed using the character property: 

dU X r(VU) X s(U^W) = L 5 Xr(VW). (10) 

d r 

In other words the junction of the surfaces in the space of 
character coefficients is represented by an ordinary prod- 
uct. 

For any 2-dimensional surface we can expand the parti- 
tion function (J5J) in characters according to © and then 
integrate all the internal plaquettes using ([9]). The re- 
sulting expression for the partition function is 



(11) 



led A 



where A = N r Nt is the area of the total surface in pla- 
quette units (number of tiling plaquettes) and dA is the 
contour enclosing the surface. 

The multi sheet n-replica partition function Z n for 2-d 
model is shown in Fig. [T] Remember that the links of 
time boundary of region A for each replica (bold lines) 
are identified and so are the links from region A that form 
the time boundary of Z n (dotted lines). Being subjected 
to the same treatment as Z the partition function of the 
glued system Z n will result in the same expression (jlip . 
but with corresponding surface area A n = nA = nN r N t 
and perimeter dA n . 

To perform the perimeter integration first we choose 
free boundary condition (b.c.) in the spatial direction. 
The invariance of the group integration (Hurwitz/Haar 
measure) allows one to manipulate the link variables in 
the perimeter integral, so that the final integration is 
performed over a single plaquette perimeter (c.f. Gross- 
Witten one plaquette integral) in both dA and dA n inte- 
grations. Specifically we absorb the space-like links that 
separate different replicas in Z n (time boundary of re- 
gion A) so that the system becomes identical to a simple 
plaquette with free spatial b.c. 

Due to periodicity in time direction the ordered con- 
tour product of gauge fields generally has the form 



(12) 



Here U n j denotes the gauge field at coordinate n in i = 
0. 1 direction, where is chosen to be the time direction. 
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We use another property of character integration 
dU sXr(U ori V h() Ul ) = j-Xr(F 1)0 )Xr(F 2 t e ) 



. (13) 

The integral over the remaining two gauge variables 
decouples and has support only for the trivial represen- 
tation Xo = 1- This leads to a simple result 



Z = F' 



(14) 



The ratio of the partition functions is unity and the en- 
tanglement entropy is zero. 

Next we consider a lattice periodic in the spatial direc- 
tion which effectively mimics infinite spatial extent. The 
perimeter integral for Z now is 

[dV [ dU Xr (UVUW) = [ dV± X r(V)xr(Vi) = i 
J J J d r d r 

where in the last part we used the character orthonor- 
mality property. The partition function becomes 



Z 



(16) 



It is easy to check that the Z n perimeter integral results 
in 



dUi. 



■duj- 
d 



Xr{U 1 )...Xr(U n )Xr{U\)...Xr(U n \) 



dr 



(17) 



Note that this expression for n = 1 correctly reproduces 
the result of perimeter inegration for Z, c.f. (TT5)) . The 
partition function ratio is 



Z n 



E r Fr A /d 2 r n - 2 

(Er F r A ) n 



1 + Er#0 C r/ d: 



nA I j2(n-l) 



(i + E 



q 4 )™ 



(18) 



The entanglement entropy then is 



dn Z n 



log(l 



r#0 



Er^O c f |QgCr/jr 

i + E 



(19) 

Note that that the series of character expansion coeffi- 
cients is vanishing (1 > c r > c s |if d s > d r ) in such a way 
that the sums in (118|) are converging even for the smallest 
surface A — 1 . One then can choose the surface area large 
enough to guarantee that (TITI)) is hnite. We observe that 
this expression is /-independent [26j and valid for I > 0. 
The entanglement entropy expression (|19[) is universal in 
the sense that it does not depend on the initial lattice 
cutoff and is dependent only on the physical dimensions 
of the system. This is due to the fact that after a number 
of iterations the coefficients are attracted to the renor- 
malization group (RG) trajectory independently of the 
starting point. It is interesting that at I = Z n factors 



into n copies of Z so that ratio Z n /Z n = 1 and Sa = 0. 
This is reminiscent of the 2-dimensional theory end-point 
phase transition at temperature T = 0. 

If the surface area is very large one can truncate the 
series to obtain a manageable expression. This is in fact 
similar to the strong coupling limit treatment. Using 
strong coupling expansion in evaluation of F r we obtain 
an approximate expression for the entanglement entropy 
(|T9|) . In general one can compute F r term by term to any 
desired order. For our purposes, however, it is enough 
to keep the first two lowest order terms, which give the 
coefficients for the trivial r — and fundamental r = 1 
representations 



dU 1 



jX* r {U). (20) 



Thus F Q = 1 and a = F t = (3/(2N 2 ) for N > 2 (note 
that characters of SU(2) group are self-conjugate and 
therefore c\ = (3/N 2 ). The entropy becomes 



S A 



log(l + (J* )*) 



2N 2 ' 



JL 
2N 2 



1 - log ( 



2N 2 



>-7-V 



(21) 



Simplifications can be also achieved in the large N 
limit. The expressions for the first two representations 
F r (integrals (O) are readily available [6j. In the Gross- 
Witten paper notation Fq = z and ci = u) 



F\ = ujz = F x 



1/A, 

1 - A/4, 



A > 2 
A < 2 



(22) 



where A = g 2 N is the 't Hooft coupling. 

Again for very large surface area it is reasonable to 
assume that the terms in this series are rapidly vanishing. 
Therefore the entanglement entropy becomes 



S A 



(23) 



Note that for the strong coupling w = 1/A = P/(2N 2 ), 
and the expression for Sa is equal to the strong coupling 
expansion derived earlier (we can interchange strong cou- 
pling and large N). It is no surprise that the entangle- 
ment entropy is sensitive to the 2D Gross- Witten phase 
transition and is different for strong and weak coupling 
phases. 



B. d > 2 gauge theory 

Next we consider the D = d + 1 dimensional theory, 
with d > 2. This is a nontrivial theory which cannot be 
solved exactly. We employ the Migdal-Kadanoff (MK) 
0, H 9 decimation procedure to solve this theory ap- 
proximately. For illustrative purposes we concentrate on 
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FIG. 2: Z n for 2 + 1 dimensional theory. 

2+1 theory, the generalization to higher dimensions is 
straightforward. 

In general for a finite temperature system one has to 
use anisotropic lattice. The time and space direction 
bond moving can be performed independently, c.f. A- 
and p- transformations for finite temperature gauge the- 
ory [HI [Hi 03 ■ To simplify the treatment we consider a 
vanishing temperature system in a symmetric box. 

The standard MK decimation procedure (A- 
transformation) moves the internal plaquettes to 
the hyper-surfaces which constitute the elementary cells 
of the resulting coarse lattice 



-S P (U) _ 



Fr 



Y,F r A d rXr (U) 



C 1 



(24) 



dUe -C b s P (U)_L xUu) . 



where the choice 6 = corresponds to Migdal, while 6=1 
to Kadanoff prescription. Here £ = \ D ~ 2 is the factor 
by which we strengthen the interaction on the resulting 
coarse lattice in order to compensate for missing internal 
plaquettes, A = A 2 is the surface of the new elementary 
plaquette in units of fine lattice plaquettes (number of 
tiling fine plaquettes), and A is the scaling factor of the 
RG transformation and is equal to the number of plaque- 
ttes (internal and from the surface) moved to the surface 
from each of D — 2 directions. 

It is known that the Kadanoff procedure results in the 
overcompensation of the strength of the coupling thus re- 
sulting in the upper bound for partition function, on the 
other hand leaving the coupling on the surface unchanged 
while dropping internal interactions £ = 1 results in the 
lower bound on the partition function [l3j 



Z(C = 1) < Z < Z{( = A 



D-2\ 



(25) 



This relation relies on translation invariance and there- 
fore does not hold for Z n , however one may expect it 
to hold approximately. As a result a generalization of 
the MK procedure which preserves the partition func- 
tion may be possible to construct. Here, however, we use 
the standard MK decimation. 



After each step of decimation iteration the partition 
function decomposes into the product of the coarse lat- 
tice partition function and the integrated out bulk part, 
which (after m step iteration) is 



|A|/A jD 



(26) 



Significant simplification can be achieved if we cary out 
decimations for Z n and Z in exactly the same manner. 
As a result of equal volumes the bulk contributions in Z n 
and Z n are identical and cancel out in their ratio at each 
step. 

We start with a symmetric D = d + 1 dimensional dec- 
imation (A-transformation J2U) in Z and Z n , see Fig. [5] 
Note that there are periodicity conditions in t-direction 
for each A part of n-rcplicas (bold links) and for the 
links of time boundary of the glued system (Z n ) belong- 
ing to A (dotted links) . The decimation should be altered 
when the lattice spacing becomes equal to I (the smallest 
scale in the problem). At this point the /-like plaquettes 
(directed along I) inside the slab of thickness I (extend- 
ing through all n replicas) have to be treated differently. 
These plaquettes can be decimated in remaining direc- 
tions, very much like time-like plaquettes in the finite 
temperature gauge theory treatment. Such transforma- 
tions are normally referred as p-transformations, for them 
the decimation prescription (|24[) is modified 



S p; i(U) 



F r 



Y.FrdrXriU) 



c 1 



(27) 



dUe 'C b S p ,(U)^ x:(u y 



We still can move plaquettes in D — 2 direction but the 
tiling is done with A plaquettes. All the other plaquettes 
are unaffected by this change and are decimated accord- 
ing to standard (A-transformation) procedure. 

In appendix A we consider a gauge theory formulated 
in a box R 3 . This system will be the building block for 
construction of expressions for Z n and Z. 

Let us assume that the imaginary surfaces that cut 
out the part for which we compute the entanglement en- 
tropy belong to x — t planes and are distance / apart in 
y direction (note that we consider 2+1). 

We begin with Z. The surface with normal along x 
consist of 3 pieces after decimation is stopped. There is 
exactly the same contribution from the surface with nor- 
mal —x. At the center there is a boundary of the slab 
c s x and two pieces which complement it, we refer to their 
joint as c x . The combined contribution is c x = c s x c x and 
should be substituted into the corresponding equation 
from the appendix. We note here that technically the 
complement to the slab (more precisely two complemen- 
tary volumes) are not symmetric, therefore the recursion 
at some point has to be switched from A to p, however we 
consider the scale R » I and therefore we can always 



5 



take R large enough so that the corresponding coeffi- 
cients are in the strong coupling limit and no transition 
from the RG flow to the infrared fixed points can occur 
in these bulks. 

The 2 surfaces with normals ±y each contribute c y . 
There is only one group of surfaces (similar to c x ) with 
normal —t (our convention) with contribution ct = c a t c 8 t . 
Therefore from (14"T1) 



C V,i) ( c x,i c x.i c y,i) djC 



After many successive steps of decimation iteration the 
only remaining degrees of freedom are defined on the sur- 
face of the system. In case of Z n we also have n — 1 i-like 
plaquettes inside the bulk (c* •). At this point we can 
move these Z-like plaquettes in Z n to the bottom sur- 
face. This decimation step has no counterpart in the de- 
nominator and therefore the bulk term (i*o) of this last 
decimation procedure does not cancel. This decimation 
step is achieved only with moving the internal plaquettes 
along the time direction onto a single surface plaquette; 
there is no integration of the tiling plaquettes for this 
procedure, therefore the new coefficients for the result- 
ing surface plaquette after Kadanoff type moving are 



XAU ] ) (29) 



there is no periodicity requirement for these links we can 
integrate them out. As a result the internal n—1 time-like 
surface terms of A have support only at the trivial repre- 
sentation and therefore do not contribute to the partition 
function. There is still, however, a contribution from the 
first replica time-like surface (bottom) of Z n . 

After simple considerations one can convince oneself 
that the surface integral in Z n is similar (in 2 + 1 di- 
mensional theory up to factor to the surface 
integral of a nN t x N% cube [13]. The side surface coeffi- 
cients are modified to account for gluing n replicas, while 
the bottom surface coefficient involves a term computed 
according to ([50]) and is cf^cf i . The partition function 
becomes 



Zn — F t S Q ■ f n — F t S Q X 

1 + y ] ,4(n-l) ( C x,i^x,i c V,i) 



(31) 



l(n-l) 

The ratio of the partition functions including the bulk 
term is 



Zn f, s 



(32) 



1 + Ei/o( C i,i C i,i C !/. 1 ) 2 ' l M' 


n-1) 




(l + J2i=io( C x,i^x,i C y,i) 2 




r 



and 



F 



I -J 



F, 



l A) 



(30) 



In order to obtain a higher dimensional expression for this 
ratio one needs to adjust accordingly the sides contribu- 
tion and the contribution from the surface integration. 

The entanglement entropy is 



Next we assume that the boundary between region A, 
which has space-like links with coordinates t = and 
t = 1/T identified and region A which has no such con- 
straint is defined in such a way that the end links (di- 
rected along x) of the cut in Z n belong to region A. Since 



S A = -F ti0 +logZ-^ (33) 



where the dot stands for X = -£-X\ , . Note that 

dn \n=l 



u y,iJ 



J2 d 



'jCtjF* 1 ^ 



(34) 



In order to manipulate these expressions we will need the The expression for the entanglement entropy can be 

following derivatives evaluated if the system flows towards the IR fixed point. 

. s ls s This is the strong coupling limit for c s t ■ therefore we can 

£ t,j = F t,j- c h F tfi (35) expand logarithms and simplify the expression. 



K,i = J M[\ + ]T dicliXiiU)} log(l + ]T di<4 iiXi {U)) 



l x ,(C/t), (36) F t ^o = 4, + £ ^clcyDi, = ct, + 0(c 2 ) 

dj M'#0 3 



6 



Ft,o = E '-<^J-''j / dU Xi (U) Xj (U) =0(c 2 ) 



The leading term in the entropy is 

Sa ~ -(c s x ,iC s xA c yil ) 2 log(c s xA c s xA c Vil ) 2 



(37) 



Note that the dependance on / is encoded in the value of 



C. Analyzing the RG flow 



Now recall that our choice of temperature makes the 



box symmetric and c: 



t,i 



. i = cf . The resulting expres- 
sion for the entanglement entropy (|33p is a very compli- 
cated function of cf. Note that this is a general fea- 
ture valid for higher dimensional theories as well. The 
dependence on I enters through the value of these co- 
efficients. Essentially I regulates the moment when A- 
transformation is switched to p-transformation, which in 
turn sets the initial value for the cf (mo) iteration under p- 
transformations thus defining where the theory will flow 
before reaching the boundary. 
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FIG. 3: Migdal decimation flow for 3+1 dimensional SU(2) 
gauge theory. Projection to cf/ 2 and cf; (f3, A) are indicated. 

Next we analyze the RG flow of SU(2) gauge theory 
for c|(m) as a function of number of iterations m un- 
der Migdal recursion (f2"T)) and depending on the starting 
point. In Fig. [3] we plot the projection of the flow (for a 
3 + 1 dimensional theory) from the infinite dimensional 
space of character coefficients onto the fundamental- 
adjoint Ci/2 — c\ plane. We consider A = 1.1 and 2 values 
and observe a significant dependence on the choice of the 
scaling factor. In what follows we will use the former 
value, since it is known to reproduces the SO (3) critical 
coupling value [13] • This value was also used to extract 
an approximately correct phase diagram for the mixed 
action fundamental-adjoint SU{2) gauge theory (l5j . 

One can clearly observe that depending on the starting 
value the flow will go to either of the two fixed points - the 



infrared trivial fixed point or non-trivial UV fixed point. 
This is a clear indication of a transition. The starting 
point for the system is set by the A- transformations and 
depends on the value of I. Generally at the starting point 
the action is a single plaquette action but with an infinite 
number of couplings for terms in all irreducible represen- 
tations. In the numerical simulation that gives Fig. [3] we 
simplify this situation by considering a starting action 
in the wilsonian (only fundamental representation) form 
on N t — 1 lattice, noting that this should not affect the 
observed picture of existence of transition. 

The lattice inverse coupling value (3* e (0.62, 0.65) 
where the transition in the flow occurs should be com- 
pared to Nt = 1 gau ge t heory finite temperature phase 
transition f3 c ~ 0.86 [Tg|. This allows one to relate the 
scale I* of the entanglement entropy transition to the fi- 
nite temperature phase transition scale l c = 1/T C . For 
this we use the standard 1-loop scaling relationship 



a(fi)Ai 



( P 



\2Nb 



bl/2bl 



exp 



4Nb 



(38) 



where b Q = 11/24tt 2 and h/2bl 
the couplings we obtain 



51/121. Substituting 



l*/l c € (1.56,1.66). 



(39) 



III. DISCUSSION OF THE RESULTS 
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FIG. 4: Migdal decimation flow for 2 + 1 dimensional 5(7(2) 
gauge theory. Projection to c\j 2 and cf; A = 1.1, /3 = 
3.0,3.1,3.2 and 4.0 (inlet). 



In this paper we studied the entanglement entropy in 
d+ 1 SU (N) gauge theory. We use the multi-replica trick 
to relate the entanglement entropy to a simple ratio of 
partition funcitons. The d = 1 theory is solved exactly. 
Free spatial b.c. lead to trivially zero entanglement en- 
tropy. Periodic spatial b.c. show non-zero universal value 
independent of the size I of the entangled region. As the 
entangled region is removed the entropy becomes zero, 



7 



showing in this manner behavior similar to the end point 
phase transition of 1 + 1 dimensional theories. 

Using MK decimation we approximately computed the 
ratio of partition functions and entanglement entropy in 
d > 2 dimensional theories. A note of caution should 
be made regarding our choice to cary out the decima- 
tion for Z n and Z in the same way. This allows us to 
significantly simplify the computational procedure. The 
non-analyticity in the RG flow observed for Z n , however, 
is also induced in Z by this choice. This should not be 
a problem if one is interested only in the location of the 
transition. 

In the case of 3 + 1 SU (2) gauge theory we demon- 
strated that there is a non-analytical change in the RG 
flow for coefficients of character expansions which define 
the entanglement entropy. We find that the length scale 
of this transition is l*/l c £ (1.56,1.66). Unfortunately 
the systematic error due to the use of the MK approx- 
imation is not easily tractable. It is interesting that in 
large N c case it was shown [2J that l* c jl c = 2. 

It is important to note that the MK procedure does 
not find a transition in the RG flow for 2 + 1 dimensional 
theories. This transition is only observed for d + 1 the- 
ories with d > 2. Most likely this is an artifact of the 
MK procedure and d = 2 theory exhibits a transition 
similar to higher dimensional theories. The MK deci- 
mation is known to miss the order of phase transition 
while correctly identifying its location. It is conceivable 
that in d = 2 theory a proper transition can be seen 
by the MK procedure as a cross-over. This observation 
is supported by the fact that we indeed observe an in- 
teresting qualitative change in the flow around (3 = 3.2, 
see Fig. |4j For values of lattice inverse coupling below 
this value the flow is directed immediately toward the IR 
critical point (monotonously decreasing series of c$, Vi), 
while for larger values of j3 the flow is directed from the 
IR fixed point for a few steps of iteration then switch- 
ing to the flow toward the IR fixed point (initial increase 
of Ci/2 followed by monotonous decrease). We illustrate 
such a scenario for weak coupling regime [3 = 4.0 in the 
inlet of Fig. [4j It is interesting that formally the MK p- 
transformation in d+ 1 theory with compact direction can 
be effectively viewed as A- transformation with an effec- 
tive RG scal ing parameter vA in 2D — 2 = 2d dimensions, 
see [13, EH El- Thus d > 2 theories are related to the 
zero temperature theories above the critical dimension 4 
(have bulk phase transition), while d = 2 is related to 4 
dimensional zero temperature theory. 

Similar results hold for SU(3) and other N c groups. 
Therefore our claim is that the transition in the entan- 
glement entropy is observed for any number of colors N c 
and the critical scale I* where transition takes place most 
likely is N c dependent and asymptotically reaches 2 as 
N c -> oo. 

We note that the finite temperature phase transition 
studies of SU(2) and SU(3) gauge theory within MK 
formalism [1(3, [ill, E3 relied on the same analysis of the 
RG flow. It is important to emphasize that the periodic 



boundary conditions in time direction do not play any 
role in such studies. One has to impose the periodicity on 
the N t = 1 system after A-transformations are switched 
to ^-transformations [l 01 ] . This will result in an effective 
lower dimensional spin system which exhibits a phase 
transition for d > 2. 

The study of the entanglement entropy effectively is 
transformed into an MK analysis of a gauge system de- 
fined with one compact direction and no periodicity im- 
posed. Possibly Monte-Carlo simulations of such systems 
can define the location of the transition more accurately. 
However, this would be still a crude approximation since 
the MK treatment results in a well defined boundary, 
which is in reality rather soft. Therefore direct numeri- 
cal computation of the entanglement entropy should be 
preferred. 

It is also interesting to relate our results to studies of 
the vortex free-energy order parameter [l7j ]. which pro- 
vides a complete characterization of the possible phases 
of gauge theory. For SU (2) it was found [l8| that when 
the transverse size of the lattice is around 0.7 fm there 
is a sharp cross-over in the vortex free energy. This 
cross-over has an obvious physical interpretation: the 
lattice size has to be large enough to accomodate suf- 
ficient spreading of the vortex flux ('fat' vortex) to enter 
the regime of exponential free-energy lowering by further 
spreading, i.e the confining or color magnetic mass-gap 
creation regime. Assuming that a/ct = A20MeV (a is the 
string tension) we get in this theory 1/T C = 0.681/to. 
Therefore the transition in the vortex free energy hap- 
pens approximately at 1/T C scale. 

Using this observation we suggest that the transition 
in the entanglement entropy happens when the size of the 
entangled region is large enough to accommodate a fat 
vortex. The difference in the geometry should account on 
small difference of the scales when such transition occurs. 
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Appendix: 2 + 1 dimensional gauge theory in a box 

In this appendix we consider a 2 + 1 dimensional 
SU(N) gauge theory formulated in a symmetric box R 3 
at temperature T = l/R, which corresponds to confined 
phase temperatures for sufficiently large R. The theory 
is formulated on a lattice with the UV cut-off a, with pe- 
riodic boundary condition in t direction and free bound- 
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FIG. 5: 2 + 1 dimensional symmetric box. 

ary condition in spatial directions. The MK decimations 
(f2"4"]l with scale factor A are performed iteratively N times 
(X N = R = R/a). At this point all degrees of freedom 
are "pushed" to the boundary and the resulting lattice 
spacing becomes equal to R, see Fig. [5] Here we are not 
interested in the bulk contribution. There are 8 indepen- 
dent gauge degrees of freedom which live on links. We 
use normals to the cube faces, which are directed outside 
to identify plaquettes. The partition function has con- 
tributions from 4 plaquettes with normals along spatial 
directions and one plaquette (due to periodicity) from 
t-direction 

f({c z };dA) ^l + ^diC-iXiidA,), (40) 

where z = ±x, ±y, t marks the plaquettes. The character 
coefficients c z -i can be obtained numerically and are the 
result of the RG flow in infinite dimensional coupling 
space. The symmetry of the box implies c x -i = c y , = c t i . 
For a general non-symmetric box one has to consider a 
series of p-transformations, resulting in an anisotropic 
lattice with all coefficients c Zj i different. 

Because of the free spatial boundary condition we can 
further integrate out three time-like links (thin lines in 
Fig. [5|. By doing this we join the surfaces according to 
©, with the resulting surface term f({c XVl i}; U'VUV') 
and c XVi i = c x { . The partition function is 

Z = [ dUdVf({c xy!i };U^VUV^f({c tji };V) 



= 1 + Cx v> 1 + r r^--' l j r '-j l Kr (41) 

where 

D$j = J WXk(V*)Xi(V)xj(V) (42) 

we recognize as the coefficients of the Clebsch-Gordan se- 
ries £> (l) x D<>') = J2k D ipW for the Kronecker product 
of irreducible representations. Using Gaunt's formula 

\Gl l f V^\R- x ) nimi V^\R) m V^\R) n37n3 dR 

JG 

= ( h \ ( h V ( ji h h V ( h h h \^\ 

where \G\ is the volume of the group space, we can ex- 
press D%j through the Wigner coefficients and 3-j 
symbols) for general group [l£j, l2fj|] 

D k = (k\(k\*(ki jY(ki j \ 

(44) 

Coefficients D r rs can be easily evaluated for SU(2) 
group, using the Clebsch-Gordan equation 

i+j 

fc=|»-i| 

Thus the integral becomes 

D h = / dVY,Xk{V) X j{V) = Hi(2i - j), (46) 

J fc=0 

where Hi(x) is the Heaviside step function (Hi(0) = 1). 
Therefore 

Z SU(2) = 1 + Y, Cx V- 1 + c xy,idjCt,j (47) 
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